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RADIAL LOAD FLOW PROBLEMS 
IN DISTRIBUTION CIRCUIT ANALYSIS 



Part A 
The Problem 



This case study covers what is, in my opinion, an 
extremely interesting application of computers to 
solve an engineering problem. The problems associ- 
ated with the distribution function of an electrical 
utility are often overshadowed by those in the 
generation and planning areas. This application 
shows that there are unique problems in distribution 
which can be solved on a computer and thereby 
relieve the engineer from relatively tedious jobs. 
This study is also interesting in that it involves 
solution of a non-linear problem. It appears that 
most textbooks only consider circuits which are 
linear in nature. However, in real life, many 
problems are non-linear and require special logic 
in order to arrive at a satisfactory solution. 

Francis A. Long 
Pennsylvania Power 
and Light Company 
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INTRODUCTION 

One day in December, 1969, Harold E. Taylor, Senior 
Engineer, Distribution Lines and Substation Department of 
Division Operations , Pennsylvania Power and Light Company , 
called Frank Long, who was Senior Programming -Engineering 
Analyst in the Methods Section of PP^L's Financial Division. 
"Frank, couldn ! t we do something to calculate single-phase 
loads correctly?" asked Taylor. 

Following this contact, a few meetings between Frank 
and Harold showed that what was needed was a more precise 
way to analyze a distribution circuit and calculate its 
voltage profiTe in order to determine the proper size and 
location of fixed and switched capacitors on PP§L's distri- 
bution circuits. A recent study made by the company had 
shown capacitor usage to be indeed an economical method 
for controlling voltage at the customer's meter in order 
to satisfy Pennsylvania Public Utility Commission regula- 
tions. 

The large, fast digital computers used by PP$L could 
handle a comprehensive distribution circuit analysis pro- 
gram which would not only help to spot capacitors at proper 
locations but aid in other ways to guide decisions neces- 
sary to control distribution voltages in an economical way. 
But Frank said, "What I had in mind was not merely a replace- 
ment for existing rough ways of calculating these things. 
We really wanted a method to reduce the manpower required 
to manage a distribution circuit. It is a non-linear cir- 
cuit analysis problem and later, when we get the system 
working we should be able to locate and size capacitors 
in the optimum manner." 

ERANK LONG , JERRY GUTH , ANJ2 ItiE METHODS SECTION OF PENNSYL - 
VANIA POWER AND LIGHT 

The Methods Section performs all scientific computer 
programming for PP§L. Equipped with the IBM 360/65 and 
360/50 computers as well as some older, "second generation," 
computers such as IBM 7074 and 1401, the Methods Section 
is composed of five engineers and eleven mathematicians, 
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three of whcr.1 are women, who program problems arising in 
PP$L f s work and send the programs into the company's com- 
puter cejite r . 

Frank Long, 1963 graduate from Northeastern University, 
has done additional studying both at PP$L and at nearby 
•Lehigh University. After five years experience as engineer 
in PP^L's System Planning Department, he came to his present 
position in the Methods Division two years ago. "With sys- 
tem engineering experience and with further work in computer 
methods, I believe I can guide programming from both the 
engineer 1 s and mathematician' s viewpoints . It's an inter- 
esting, creative job, and I enjoy the challenge," Frank 
remarked . 

After the initial discussions with Taylor, Frank roughly 
estimated the job would take a certain number of man-hours 
and selected Gerald Guth, a scientific programmer, to work 
with him. Jerry graduated from Lehigh University in 1966 
with a degree in engineering physics and later a master's 
degree in math. At the time of this distribution circuit 
assignment, Jerry was working on a master's degree in phys- 
ics from Drexel Institute of Technology under PP^L's tui- 
tion repayment plan. Frank said, "I asked Jerry to work 
with me on this problem because I felt he was capable and 
had the time available." Jerry came to the Methods Section 
in 1968, after working in the Atomic Power Section for two 
years since joining PP§L. 

"The job," Frank said, "was actually not of a top pri- 
ority nature. We'd had something that worked for some time 
which we could live with--an approximate way of treating 
the case. However, we all felt that in the end significant 
improvements would result." 

Frank further explained that Jerry and he had allot- 
ted time to work on this job steadily unless urgent jobs 
came in- -and there had been several such instances. "For 
example," Frank pointed out, "Jerry once had to drop every- 
thing in order to check and correct an ASME program which 
calculated enthalpy from pressures and temperatures. This 
work was done on a hurry-up basis since our Brunner Island 
No. 3 generation plant badly needed the correct data." 



PENNSY LVANIA POWER AND LIGHT COMPANY 

The company serves a 10,000 square mile area in Central- 
Eastern Pennsylvania. It is an investor-owned public util- 
ity wi'h 98,450 shareowners; its headquarters is in Allen- 
town, Pennsylvania. PP£L is part of a large power pool 
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known as the Pennsylvania -New Jersey-Maryland Interconnec- 
tion, composed of 13 companies serving an area of 48,700 
square miles. PP$L produces electric power for residential, 
rural, commercial, and industrial users from its five steam, 
three hydro- electric, and nine combustion-turbine powered 
generating stations and, in time of heavy demand--in winter- 
time, at present, due to an increase in number of electrically- 
heated homes--can also purchase power from neighboring util- 
ities . 

In 1959, peak demand was 1292 megawatts; in 1969, 2702 
megawatts- -more than double. The company management expects 
the demand to continue increasing, and in fact, lists the 
following as major problem areas: increased taxes, infla- 
tion because of higher costs for construction and operation, 
increased demand for electricity, environmental protection, 
and higher standards for electric service reliability. 

Figure 1A of Appendix A shows the main PP^L system 
in 1974. Figure 2A gives further data on PP$L. Figure 
3A is a simplified organization chart with the positions 
of the principals in this case shown. 



LEGAL AND PRACTICAL ASPECTS 0£ THE PROBLEM 

The Pennsylvania Public Utility Commission sets oper- 
ating and finance standards for all investor-owned electric 
utilities in the state. The regulations specify that volt- 
age (60H ) at the customer's meter shall not exceed the 
z 

nominal standard (120 volts for PP$L) by ± 5%, or between 
114 and 126 volts, and, further, that the total variation 
between sunset and 11:00 p.m. shall not exceed 8%. Thus, 
one objective is to comply with the regulation, but also, 
as Frank commented, to do so in an economical, reliable 
manner . 



UNDESIRABILITY OF VOLTAGE VARIATION 

Frank further pointed out a study done by another PP$L 
engineer on the effects of voltage variation on customer 
utilization equipment. This study gave figures to support 
the common-sense conclusion that the best combination of 
reliability, lifetime, efficiency, and economy is achieved 
if electrical equipment is operated at the design voltage. 

For example, induction motors overheat or may not start 
at all, resistance heating devices take longer to heat up 
on the one hand or shorten in life on the other, and tonal 
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quality and volume, brilliance, and focus of black-white 
and color TV sets may be affected. Lumen outputs and the- 
oretical life of incandescent lamps are drastically altered 
by voltage changes. Figures 4A, 5A, and 6A of Appendix A 
suggest tiic importance of maintaining proper voltage. 



THE DISTRIBUTION CIRCUIT MANAGEMENT PROBLEM 

The distribution circuit, as explained by Frank , con- 
sists of the power lines from the substation to the cus- 
tomer's location. t! It usually begins with a 66-12KY step- 
down transformer and is made up of a main three-phase line 
out to a certain point, with branches of three-phase or 
single-phase taps from p lace - to -place , ,f Frank said. (Figure 
7A of Appendix A shows an actual lay-out of a circuit.) 
Wires are strung on poles according to standard techniques, 
two of which are shown for 12KV lines in Exhibit 1A of 
Appendix A, page 9A. 

Tnc input data for a distribution analysis consists 
of the KW and KVAR--+ for inductive, - for capacitive- - loads 
on each phase, the phase connection data, distances between 
nodes, wire sizes, and voltage at the source, or bus. Us- 
ing the wire data and construction, or pole, type, the resis- 
tance and impedances per mile may be calculated (reference 
1 in the bibliography). 

Figure 8A of Appendix A gives conductor data typically 
used and Exhibit 2A gives a printout of circuit data for 
an actual circuit on the PPqL system. With this data, a 
voltage profile for peak and light load conditions, with 
and without fixed or switched capacitors on the system, 
can be plotted to be sure that voltage criteria are main- 
tained at all times. 

Capacitor banks, connected at various points, are used 
to control voltage levels by providing negative KVAR's to 
loads requiring positive KVAR's. Frank explained, M It T s 
the power factor correction problem, actually. But the 
circuit equations are non-linear since we assume a constant 
KVA load--a little bit different from the problems you find 
in textbooks . t! 

Primary distribution systems should be designed and 
operated so that transformers can be economically loaded 
and still provide a customer voltage variation within allow- 
able limits. On the average, with some allowance for over- 
loading the distribution transformer, the voltage drops may 
be apportioned as follows: 
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Peak Load 



Light Load 



Drop from 66KV to 12KV 
bus 



.5% 



Drop in 12KV 1 ine 



1% 



.9% 



Drop in distribution 
trans former 



3% 



1 .3% 



Drop in secondaries 
and services 



3% 



1.3% 



Total 



9% 



4.0% 



According to Frank, the preceeding apportionment does 
not make efficient use of the system facilities, particu- 
larly in the 12KV circuit unless voltage compensation is 
applied. For example, if limited to a two percent voltage 
drop on the 12KV circuit, a 10-mile, 3-phase line of No. 2/0 
copper conductor, uniformly distributed load among the 3- 
phases, would be limited to less than 10% of thermal capa- 
bility . 

On the PP$L system the voltage at any given location 
can be established by the use of one or more of the follow- 
ing alternate methods: 

For general level control--a rough adjustment-- 

1. Taps on the 66-12KV transformers 

2. Fixed capacitors on the distribution circuit 

3. Taps on the distribution transformers 
For voltage compensation- -defined as 



1. Variation of the 66KV voltage by varying the 
generator field strength, and tap changing 
under load of the 220-66KV transformer in 
the substations 

2. Switched capacitors on the 12KV distribution 
system 



11 a.m. voltage - 3 a.m. voltage 
3 a.m. voltage 



x 100-- 



3. 



Use of voltage regulators. 
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In 13C2 PPfiL concluded from a broad study (see refer- 
ence 3 in the bibliography) that the optimal method for 
controlling distribution voltage was by the application 
of fixed and swi tched capacitors -- fixed capacitors would 
be connected at a given location and left for a longer pe- 
riod of time; for short-term needs, other switched capaci- 
tors would automatically switch on during peak load times 
and later switch off the lines. Application of this pat- 
tern is intended, Frank said, to maintain a primary distri- 
bution voltage within the limits of 120-126 volts during 
peak load periods and 119-123 volts during light load pe- 
riods. These voltages are not the actual primary voltages 
but the primary voltage referred to secondary. 

Frank said that the purchase of variable tap trans- 
formers was discontinued in the spring of 1962 because a 
change in the manufacturer's transformer standards added 
a premium of $20.00 for taps to the purchase price of a 
10 KYA transformer and $37.00 to a 25 KVA transformer, 
for the PP'ul Company, this would have resulted in a total 
increase of $132,000 for a year's supply of transformers. 

Also, the use of voltage regulators on the distribu- 
tion circuits was undesirable from both the consideration 
of initial capital outlay and later high maintenance costs. 

Frank pointed out that the use of capacitors, properly 
sized and located, not only promised better economy but 
also would improve the over-all system power factor and 
thereby release thermal capability on the entire system 
from the generators to the distribution transformers at 
significant savings . 

For example, the 1962 study on use of capacitors stated 
tnat using $3 million to supply capacitors to the Plh'il. sys- 
tem would raise the total system power factor at peak load 
from 80 to 90i, possibly resulting in decreased transmis- 
sion and distribution losses saving, under optimum condi- 
tions, $42 million/year! Actually this large saving esti- 
mate was not realistic, but the re' was. no doubt at PPqL of 
the economic benefits. 



CALCULATION OF V OLTAGE PROFILES 

Since voltage is the item of interest along the dis- 
tribution circuit, plots of voltage vs. distance from the 
66 KV bus were calculated by PP§L distribution engineers, 
Frank related, first by hand and later by computerizing 
the hand method. Using the data from the simple distribu- 
tion circuit and the phasor relations shown in Figure 10A 
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of Appendix A, the voltages along the circuit were calculated 
fay Frank and Jerry using the hand, or approximate, method. 
Results are shown in Figure 11A for the circuit when the 
capacitors are connected. These numbers and the resulting 
voltage profile were to be compared with the exact method 
to be programmed. 

Examples of actual voltage profiles calculated by a 
PP$L engineer are shown in Exhibits 3A and 4A of Appendix A, 
pages 14A and ISA. 

Frank explained that, in brief, the steps commonly 
used in checking a circuit for proper voltage were: 

1. Calculation of a light load voltage profile with 
no capacitors. If the light load 12KV bus volt- 
age--120 volts used as a base--was below 120 volts, 
the 66-12ICV transformer tap change was calculated 
to bring the bus to 122 volts. 

2. Then, after attaching all fixed capacitors, cal- 
culation of a new profile for the three-phase main 
line. If the lowest voltage was less than 118 
volts, the fixed capacitors are determined to raise 
the lowest point to 119 volts. 

3. Repetition of step 2 for peak load conditions. 
If the louest voltage were found to be less than 
119 volts, additional switched capacitors were 
added beyond the lowest point toward the end of 
the line to raise the minimum voltage to 120-122 
volts . 

4. Repetition of steps 2-3 for single-phase taps. 

Other considerations entered into the problem too, 
but Frank emphasized that the key was an ability to analyze 
the circuit . 

The V = PR + QX relation was admittedly approximate 
and for single-phase loads on a single-phase circuit an 
approximation led to use of a multiplier of 4,5 to compute 
the voltage drop. Also source voltages for the single-phase 
taps were those obtained from the averages of the three- 
phase circuit voltages at the point of tap. However, with- 
out the computer to solve a large distribution circuit ex- 
actly, the hand calculation work was completely out of reason, 
so the approximate method was obviously the only solution 
available in prior days, 

Frank said, "Traditionally, solution of unbalanced 
three-phase circuit problems have been done by the method 



8 



ECL 194A 



of symmetrical components. This is impractical by hand 

here, so approximate methods had to be used. Also, the 
whole thing depends on information about the circuit- - loads , 

wire sizes, and so f or th - -which , of course, are known only 

to a certain accuracy . However , there appeared to be cases 

in actual practice, according to Harold Taylor, where a 

more accurate way, made possible by the computer, could 
be used p r o f i t ab 1 y . 1 1 

With these considerations in mind, Frank and Jerry 
began to study the problem in detail. They charged their 
own and computer time to an assigned account number and 
eventually Taylor's department will be charged with the 
development cost. 

What would your approach to this problem be? 
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GENERAL EFFECT OF VOLTAGE VARIATION 
ON INDUCTION MOTOR CittPACTERIGTICS 
CONSTAliT TORQUE LOAD 



NOTE: THESE EFFECTS WILL VARY SOMEWHAT FOR SPECIFIC EATINGS 
Figure 4A 
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SOURCE: STANDARD HANDBOOK FOR 
ELECTRICAL ENGINEERS 
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EXHIBIT 1A 



ECL 194A 



CO 




For derail of pin +i^pa- 
tnsula^ or assembly sec 
G - 13 -13 




<4 



14- 



f / 

V5B 



Jo 



/ 
"1; " 



METHOD 1 



METHOD ?. 



note: 

WHEN KtqUIRED SY OVERHANGING TREES, ■• 3/8" H.S. G'- l/2' U.S. CVERWEAD SHIELD WIRE, 
SPECIFIED BY CKCIMEEA, SHALL Be I NSTALLCO AT TOP OF PCLk PCR 6-23-36. 



For B/M set 







CV-OSE -SPKCeO PRIMARY COUDUCTOR 
PLAsSTlC 8RA.CKUT - THREE 
STRAIGHT LINE CONSTRUCTION 
12 KV AUD BULOW 
AWGLf-S 0°TO 3° 


ALLEMTOV 


$ LGT. C< 
VH, PA. 


RCV.-SKOMSOR 


/!«■•(- 


onis.AVPR.ovr:: 


• • 


1 







9A 



?.F. 



COFFaiR CQiw-JCiOR VOLTAGE DROP TABLE 

ECL 194A 

KVA MILES FOR lj VOLTAGE DROP - BALANCED 3 PHASE CIRCUIT 
12.47 KV EQUIV. SPACINc/^-8^ ) CM Q 

I/O 



100 663 1030 1610 2040 2560 

95 630 926 1350 1630 1940 

90 630 910 1300 1540 1800 

85 642 910 1280 1500 1720 

80 658 926 1270 1470 1680 

70 700 962 1280 1460 i64o 

60 752 1010 1310 1470 1630 

50 827 1080 1360 1510 1660 

to 925 1170 1440 1570 1700 

35 980 1230 1480 1610 1730 

30 1050 1290 1530 1650 1760 
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g IMPEDANCE PER MILE ON 10 MVA BASE AT 12.47 KV 

R 15 9-8 6.2 4.9 3.9 q 1 ? c , Q 

x 5-2 5.0 4.9 4.8 4.7 hS tl h9 



2 16 11 7.9 6:9 6.1 £5 1:1 ta 

KVA OF CAP AC ITORS FOR ijo RISE PER MILE OF LINE 
1920 2000 2040 2080 2130 2175 2270 2270 
KVA MILES FOR lj VOLTAGE DROP - 4„l6 KV EQUIV. SPACING 2' -10 " 



100 72 114 

95 69 104 



453 572 
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Figure 8A 
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Phasor Diagram 



E 2 = (V + AV) 2 + (.6V) 2 

= (V + RI cos <{> + XI sin <f>) 2 



( v~ ~ } 

Assuming 6V <<V + AV 

e = v + PR * 

e - v = av = HLfS* 

Also assume V Z E Z 1 p.u. 

AV = PR + QX i£ P, R, Q, X are p.u. values. 



Figure 10A 
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Approximate Hand Calculations for Figure 7A Circuit 



(KW and KVAR) 

12,000 or 120 v = 1 p.u. 
All R's and X's on 10,000 KVA 
base 

Let Bus voltage = 126 volts 
= 1.05 p.u. 



= (. 2120) (.062) + (.0025) (.092) 
= .0132 + .00023 = .0134 

E 2 = 1.05 - .0134 = 1.0366 (124 volts) 

A 23 = (.0240 x .062 + .008 x .049) x 4.5 

= (.00149 + .000392)4.5 = (.00189)4.5 = .0085 

E 3 = 1.0366 - .0085 = 1.028 (123 volts) 



P 12 = 


2120 


+ j 25 


P 23 = 


240 


+ j 80 


P 24 


1430 


- j280 


P 45 = 


750 


- j240 


P 46 = 


500 


- jioo 


AV 12 


■ P 12 


R 12 * Q 



AV 24 = (.1430 x .093 - .0280 x .138) 
.0133 - .0037 « .0096 

E 4 = 1.0366 - .0096 = 1.027 (123 volts) 



AV 45 = (.0750 x .039 - .0240 x .047)4.5 

= (.00292 - .00113) x 4.5 = (.00179) x 4.5 = .00805 
E 5 = 1.027 - .00805 = 1.018 (122 volts) 

AV,, = (.0500 x .300 - .0100 x .047) x 4.5 
4o 

= (.0150 - .00047)4.5 = (.0145) x 4.5 = .0652 

E, = 1.027 - .0652 = 0.962 (115 volts) 

o 

Figure 11A 
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Engineering Case Library 



RADIAL LOAD FLOW PROBLEMS 
IN DISTRIBUTION CIRCUIT ANALYSIS 

Part B 

The Engineering Solution 



ECL 194B 



CONSIDERATION QL ALTERNATIVE SOLUTIONS AND APPROACHE S 

"I felt the most logical thing to do was to go to the 
literature, and I found most of the really useful stuff 
in the IEEE Transactions on Power Apparatus and Systems , " 
Frank relatecT 

From the beginning, both Frank and Jerry realized the 
problem was a non-linear circuit analysis problem, as well 
as one involving an unbalanced three-phase system. The 
PP§L policy to control voltages by capacitors was firmly 
established as economically secure; therefore Frank felt 
no inclination to do other than search the literature for 
radial load flow calculations on unbalanced three-phase 
circuits . 

"There wasn't too much in the literature on distribu- 
tion circuits. However, I found an article by a group at 
Baltimore Gas and Electric which looked good, M Frank added. 
(See reference 2 in the bibliography.) Frank obtained and 
studied this article and then gave it to Jerry to study 
and begin writing a suitable computer program. 

u We might have obtained the program they used, but 
much time had passed, computer languages had changed--we 
program mostly in PL/1 here at PP§L--and it is always dif- 
ficult to decipher someone else's program, Frank explained 
He added, "Actually, the program logic is not so difficult. 

According to Frank, the basic computation scheme is 

this : 

1. Load admittances are computed using load data and 
assuming a voltage at every load equal to the 12-K 
bus voltage . 

2. These load Y's and line Z T s being used, the cir- 
cuit problem is solved once using as a source volt 
age the bus voltage, collapsing nodes and wires 

to form a new node. Nodes farther out from the 
source bus are collapsed step-by-step back toward 
that bus. Finally, the circuit is expanded out 
again to find new values for load voltages . 
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3 The load voltages which result are then used along 
with the constant load KW and KVAR to r e-compute 
a new load admittance. Then step 2 is repeated. 

4 This iterative process, alternately collapsing 
and expanding the circuit, is continued until the 
load voltages calculated in an iteration change 
less than some specified tolerance away from the 
immediately prior iteration. 

The problem is one in steady-state AC circuit analy- 
sis but with inputs of constant load power, line imped- 
ances, and bus voltage, the equations are non-linear and 
various methods could be used, among them this iterative 
technique . 

Frank further looked over books on circuit analysis 
for review purposes. Standard load flow techniques such 
as the Gauss -Seiuel or Newton -Raphson methods were discarded 
as these do not exploit the presence of a radial network 
so that a special radial load flow program could run faster. 

"I gave this paper by the Baltimore group to Jerry, 
who I felt could do a good job of programming the problem, 
Frank said. 



WRITING THE CO MPUTER PROGRAM 

Jerry related, "Well, the first thing I did was to 
read the paper by the Baltimore Gas and Electric group. 
I had to get a complete understanding of the equations used. 
Not everything in the way of needed background was in that 
paper so I had to check some of the paper's references, 
too. You have to be correct in all details if you re go- 
ing to put the problem on a computer and get correct an- 
swers. So I had to spend some time studying basic circuit 
theory ." 

A production program such as this one was to be would 
not be written by engineers in other parts of the company. 
Company policy was to assign the problem to a professional 
programmer, who could write a more efficient, more flexible 
program, it was thought. 

Jerry said that since he was not an electrical engi- 
neer and this was really his first electrical engineering 
problem he had to spend some time building up his own back- 
ground. Understanding the subject of line impedance cal- 
culations and the effect of multi -grounded neutrals with 
earth returns was the hardest part. Momentary confusion 
was caused by mis -numbered figures in the Transactions paper, 
too . 
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The CPS (Conversational Programming System) available 
at PP§L was used by Jerry to begin programming after the 
general theory of the equations to be used was understood. 
The CPS enabled Jerry to write the initial program state- 
ments in the direct interactive mode, by sitting at a key- 
board and getting instant acceptance or rejection of pro- 
gram statements typed in. 

"I followed my usual practice of breaking the problem 
down into what seem to be major parts. Then I consider 
programming each part in turn and assure myself it is tech- 
nically and logically correct, has available all its nec- 
essary input data, and is really computable. I do this 
before I try to put everything together in one package, M 
Jerry said. 

According to Jerry, the reason for this procedure is 
not only that these sub-units are sometimes best put into 
subroutines, out also he has found by his own experience, 
and that of others, that there is a danger in writing large 
programs of spending too much time on one or several parts 
and then finding out much later that a latter segment can- 
not be written or can be written only by radically changing 
all prior work . 

"Actually," Jerry added, "I didn't even write a flow- 
chart for this problem, except maybe mentally. It isn't 
that big a problem." (The brief flow chart which Jerry 
later came up with and his own explanation of the procedure 
of computation are in Exhibit IB of Appendix B , pages IB- 
4B. ) 

Jerry chose to use nodal equations for his basis, and 
let the grounded neutral terminal at the 12KV,Y— connected 
secondary be the over-all distribution circuit reference 
point. All voltages are calculated with respect to that 
point . 

Although an iterative procedure has within it the pos- 
sibility of converging to an invalid solution, Jerry felt 
it was unlikely in this case since starting values of load 
voltages are really so close, numerically speaking, to the 
true voltages . 



PROGRAM DE TAILS AND LANGUAGE 

The program was written in PL/1 language. Frank ex- 
plained this was the current language for all programs to 
be used internally by PPqL . "PL/1 is definitely user- 
oriented, and that's why we use it, but actually if speed 



4 



ECL I94B 



is all -important, Fortran is better," Jerry said. He added 
that Fortran IV is the language used if programs are to 
be widely circulated. (See Exhibit 2B beginning on page 
8B in Appendix B for a listing in PL/1.) 

Jerry also had to develop a node-numbering scheme so 
that a node closer to the source would not be "collapsed" 
too early. An end-order binary tree data structure traverse 
scheme followed by a pre-order traverse was the result. 

Although variables and computations can be declared 
and done in the TYPE COMPLEX mode, this was not done here. 
Frank explained that breaking down the work into real and 
imaginary parts and writing the program that way saves com- 
puter time in the long run. "The complex routines aren't 
'in-line' functions so they have to be called up every time 
they're needed. Programming effort is saved, but computer 
time is lost," Frank pointed out, "and since the program 
is to be used over and over again, we feel it better to 
do things this way." Double precision was not used since 
the input data was not actually precise enough, time of 
computation would go up unnecessarily, and round-off error 
was not expected to be a great problem. 

The major numerical method used is Gaussian elimina- 
tion to solve a system of eight simultaneous linear equa- 
tions, as Jerry points out in Appendix B, Exhibit IB. 

The conversion of the admittances at a node and the 
line admittances to that node into driving point admittances 
at the adjacent node nearest the source bus are done by 
means of an algorithm mentioned by Jerry and included in 
Supplement IB, page 5B. The iterative technique is also 
stated there. 



PROGRAM CHECKING AND RECONSIDERATION 

Jerry said, "There were, of course, many minor alter- 
ations as I went through the task. After it was done, I 
felt the basic computational part could be checked with 
a problem I could solve easily by hand. I tried to do this 
and got some very peculiar answers." 

This trouble caused much "head-scratching." He could 
find no mistakes in his logic. Repeated trials on problems 
with known answers gave negative results ! 

The matrix equation on page 2B of Appendix B is YE = 1 
where Y is a 4x4 matrix, E and I both 4-element column vec- 
tors and all elements are complex numbers. Jerry had written 
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(Y R + j Yj).(E R + j Ej) - I R + j Ij - 
and multiplying, separating real and imaginary parts, 

Y R E R " Y I E I = I R d) 

Y R C I + Y I E R = h (2) 

Multiplying (1) by Y R , (2) by Yj, and adding, Jerry had 
written 

< Y R 2 + Y I^ E R = VR + Y I I I 
and next, (1) by Yj, (2) by Y R , and subtracting, 

(Y R 2 + Y : 2 ) Ej = Y R Ij - YjI R . 

Thus the real and imaginary parts of the voltage phasor 
could be computed separately, using 4x4 matrices, saving 
computer time. 

M I had originally scribbled all this down on a yellow, 
coffee-stained piece of paper and was quite pleased to find 
a way to uncouple the calculations for the real and imaginary 
parts of the voltages, " Jerry related. 

"However," Jerry said, "I finally began to doubt my 
math rather than the program logic and saw I had done this 
separation by assuming that matrix multiplication was com- 
mutative! After seeing this, I corrected the program and 
answers came out all right. It's easy to get locked-in 
mentally on some things." 

A main check on the correctness of the computations 
now is a check on the power balances, real (P) and reactive 
(Q) . All P's and Q f s are computed and found to check with 
that supplied to the network. 



FUTURE COM PLETION OF IH£ PROGRAM AND USE, IN PRACTICE 

Jerry pointed out that the program is not totally com- 
pleted yet. Remaining are input-output procedures. Frank 
also added, "We want to get to the point where off-nominal 
transformers can be included. And then, of course, now 
that the circuits problem is solved, we want to study the 
problem of placing capacitors at optimum locations, and 
this is a problem in linear programming which remains to 
be done . " 
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In actual use, data cards with distribution circuit 
loads and other essential data will be maintained and when 
analyses are required by distribution engineers, the data 
can be simply fed into the computer and the program called. 
Frank said he hoped the various distribution circuit data 
could eventually be stored on disks to eliminate use of 
card input. 

"It's true that in many instances the former way of 
calculation is accurate enough, and there's always the prob- 
lem of doubt on the accuracy of basic input data. But with 
this program, we've eliminated one more possible area of 
uncertainty from the problem of managing a distribution 
circuit, M Frank said. 

To back up this point, Frank and Jerry have computed 
with the new program the circuit voltages for the simple 
6-node system, given in Figure 7A, Appendix A. Table IB 
of Appendix B, page 7B, presents their results compared 
with voltages given by the former approximate method. Jerry 
pointed out, "This system isn't too realistic, but it does 
demonstrate the ideas. The hand calculations using approx- 
imate methods give very erroneous results due to the increase 
in neutral voltage caused by the heavily unbalanced phases." 

The computer used was an IBM 360/65, which solved the 
6-node problem in 1.47 seconds after six iterations . The 
typical distribution circuit might have 100 nodes. 



ECL 19 4B 



APPENDIX B 



ECL 194B 



EXHIBIT IB 

(By Jerry Guth, Pennsylvania Power and Light Co.) 



CAL. LINE 
ADMITTANCES 
FROM INPUT 
DATA 



INITIALIZE 
ALL NODE 
VOLTAGES 
WITH SUB 
VOLTAGES 



COMPUTE LOAD ADMITTANCES 

FROM NODE VOLTAGES, 
KVA'S, AND POWER FACTORS 



INODE = N~] 



FORM DPN ADMITTANCES AT 
NODE = FROM NODE (INODE) 
BY STAR TO MESH 
TRANSFORMATIONS 



ADD DPN TO LOAD 
ADMITTANCES TO 
OBTAIN LUMPED LOAD 
AT NODE 



INODE = 


• INODE 




-1 




D 



COMPUTE NODAL VOLTAGES 
AT NODE = INODE FROM 
VOLTAGES AT FROM NODE (INODE) 
BY SOLVING NODE EQS . 



FIND A VOLTS 
FROM PREVIOUS 

ITERATION, 
SAVE LARGEST 
A IN 6 



INODE = INODE 
+ 1 




IB 



SOLUTION COMPLETE 
COMPUTE LINE LOSSES, 
CURRENTS, ETC. 
PRINT OUTPUT. 



From Node (I NODE) 
E 
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I -NODE 




Equivalent Circuit Node to Node 



Node 



equations for solving E^, E2, E^, E^ are 



VVV Y P 



•Y. 



Y 6 +Y 3 +Y 5 +Y p 



-Y 



-Y. 



-Y, 



Y 4 +Y 2 +Y 5 +Y p 



-Y, 



1 3 2 n 



E„ 



1 '3 '2 

As a matrix equation write 

YE = I 



The two four dimensional matrix equations can be written 
as one eight dimensional equation thus: 



ErY 

5 p 



E , Y 
6 p 



E 7 Y P 



" Y x • 


• V 


• 


~ E x 




1 

X 


Y 

y 


Y 




E 

y_ 




I 

_ y_ 



It is apparent that this system is no longer symmetric 
as the original complex one was. This rules out the use 
of special methods such as Cholesky's method in solving 
the system. In fact Gaussian elimination is as good as 
any, and it is the method used. 
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Computations 
Block A: 

Given KVA loads and corresponding power factors, and 
node voltages find the admittance . Let 

AE = voltage across load 
Y = load admittance 
I = current thru load 

P = KW + j KVAR = AE'I* = complex power 

KVA = |P| 



so 



KW = KVA* power factor 
KVAR = KVA* 1 - (power factor) 2 

I - AEY 



so 



P 
Y 



AE| 2 Y 

KW . . KVAR 
7 " J 



AE 



AE 



If power factor is negative take negative sign, else positive. 
Inductive load O negative power factor. 



Blocks B and C: 

Given two-node three-phase network find the driving 
point impedances at the first node looking toward the second, 
method used is to eliminate all 4 nodes of the second 'node' 
by 4 star to mesh transformations, one at a time. The gen- 
eral transformation is 




n 

E Y 
k=l 1 
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This is conveniently done for our case of eight to 
four node transformation, by the algorithm given in the 
original paper. 



Block D: 

Given the circuit of page B, with known admittances 
and voltages E 5 , Eg, E ? and E g , find voltages E^ E 2 , E 3 

and E^. 

The method used is to solve the system of node equa- 
tions of page 2B. This gives the voltages immediately, 
rather than currents as in solving loop equations. 

The equations are complex. They could be solved using 
complex arithmetic by machine, but this is very inefficient 
It is best to reduce the problem to real variables. This 
can be done by enlarging the system of four simultaneous 
complex linear equations to eight simultaneous real linear 
ones in the following manner. 

Write the complex matrix equations as 

YE = I YE - YE =1 

xx y y x 



YE + Y E =1 

x y y x y 



where 



Block E: 



Y = Y x + j Y y etc 



The voltages obtained from step D are compared to their 
values of the previous iteration. Define for the i** 1 node 



A . -max 
1 





1 




l(H»- 


1 


e^/eJ 


(ES" 


1 

• 


e^/e^ 



\ 



where E^ is the value of 

voltage E ? obtained in the 

th 7. 
n i teration . 



When A^ has been found for all nodes, then find 6 = max 



If 6< tolerance the interation is ended, 
fined by user thru input. 



Tolerance is de 
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SUPPLEMENT IB 



NETWORK DRIVING POINT ADMITTANCE AND. ITERATION METHOD 



Node 



Node 




Form array w(n,m) 



n. m 


1 


2 


3 


4 


5 


6 


7 


n ^ 
1 


y 4 » i 


y 5 ,i 




0 


y 4 »j 


y 5 ,j 


y 3 ,j 


2 


ye* 1 


y 2 »i 


0 


0 


y 6 ,j 


y 2 ,j 


0 


3 




0 


0 


0 


y^j 


0 


0 


4 




V 1 


V 1 




0 


0 


0 



Now for: 



r = 1,4 
s = 2,4 
t = s to 4 
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find 



(W 



t-s + 1 , r 



.) (W 




W 



s-l,t-s+l+r 




+ 



1, t-s+l+ r 



s 



and let 




3,5 



2,6 



1,7 




1,5 



1,6 



2,5 



Ne two r k DPN Algorithm 

These Y's are added to the Y's already existing at 
node j due to loads on node j or due to some other node 
which has been M collapsed M onto node j in the same manner. 

Repeated application of this method eventually works 
back to the first node away from the source voltage. Solu- 
tion of nodal equations here finds node voltages. Then 
working outward away from this node by repeated solution 
of nodal equations solves for voltages at each node until 
the end of the branches are reached. 

See Reference 1 of the bibliography for computation 
of phase and neutral equivalent impedances. The earth re- 
turn is not considered as zero impedance. Thus all grounded 
neutral points are not at the same potential. 

Y - 0 for an open circuit (Z = °°) • This enters zero's 
into the nodal equation matrices . 
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TABLE IB 



Results for 6-Node Circuit 



Hand (approximate) calculations : 
with capacitors added 



120 v = 1 p.u. ; 126 v 
assumed at A,B,C phases 
to-neutral at source, 
120° out of phase . 



Node 



Phase 
A 
B 
C 



2 

1.037 
1.037 
1.037 



3 

1.028 



4 

.027 
.027 
.027 



1.018 



Case I (computer results - no capacitors) 



Phase 
A 
B 
C 



2 

1.044 
1.036 
1.053 



3 

1.040 



4 

1.043 
1.017 
1.056 



1.011 



Case II (computer results - 
Phase 2 3 
A 1.047 1.044 

B 1.039 

C 1.050 



with capacitors) 
4 5 

1.052 

1.026 1.022 
1.049 



6 

0.962 



6 

1.011 



6 

1.025 
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EXHIBIT 2B 



S7*T LEVFl "F-ST 

} PL ■■:!*: J * - ,r OPTIONS* ? 1A IN) ; 

? 1 !)FC L A p E 



/*- 



STnI^LF DECISION FLOaIINS ^1^1 ST/.MC VA'MAriL 



( V 1 , V2 , V 3 , V4 , PHASE t LOO 7 » 
OELTA.E ( ftt2 I ,WURK<4 , 7,2 ) ,CUR ( ) , v f *. f /, , 2 \ » 7. P> t d ) , 
VOLTS (3^,4, 2) ,KVA(3>;< ,6) ,i;PMY ( 5 , 6 , 2 ) , P HAS FY ( 3'. t2 >, 
POwfcR_FACTnR(3U ,6>,NEUTY<3 ' f 2 ) , *> iv )f> ( 4 , <t ) , 
PF, SINPFtESUMtEXtEY ,KVA/, , 
SU-'X, SUMY t SUMSOf 

W0"*K1X,W0RK1Y , WORK2X , n'o? K2 Y , *?QPK 1 ? X , Wu" K 1 2 v , PkODX , PPCRY t 
Yl ,Y2,Y3,Y4,Y5,Y6,YP,YN, 

VOLTX,VOLTY,YPX, YPY, SI Iz ,°IVOT ,. M ULT,MAX, 
GKM f 01 V t TOLERANCE t DI FFERENCE t VIJ » 
V( S> ,DPN<6> 

> STATIC FLOAT, 

/* DOUBLE PRECISION */ 

SUM STATIC FL0AT<16), 

/* OTHER */ 

TL ST STA1 1 C F IXFD B I N ( 15 ) # 



/*- 



ITERATE: 



HAL 1 40*0 SINARY INDiClKS ANO/uR COUNTERS 

*/ 

( t , J, IY, iNODEt JNOOE tLAST.M-JOc , K » K t i- i , S , S Ml , S MP 1 , n , 
Pi (8) ,PT!,SAVE_PT,PtVOT_lNDEX t 
FR0<4_NQDE< 300) » 
T , T T , I ROW » I COL , I IROW , I PROD » t t , t t I 

> STATIC f IXCIi UINl IM, 

/* OVFKLAYS */ 

L"ADY<30e t 3i2> FLOAT PfcTlNFU UP NY ; 



TFST = C; /* SFT TEST SWITCH F OR CCNVERGFNCE CRITERION */ 



THIS ROUTINE WILL CALCULATE LUAU ADMITTANCES ASSUMING CONSTANT 
ML LuADS A T -ilVEN PJWt* FACTJR.S Ubl^ NOOE VGLTAGFS COMPUTED 
IM P*- lVI'JUS I T-RATi 0N 3 
*/ 

O'J T.^J'lOt-2 TO LAST_NU*>r. ; 
/*^«ft*#-- S rt C^LC'JIATL VOLTAGE M."-MSS Y F& \" • V)f Vr-LTAG^S #/ 



s 


1 


i 


no 1=1 to 2 ; 


6 


1 


7 


Vl^VULTSC I NOOE t 1,1); 


7 


1 


7 


V2=V0LTS< I NOOE, 2, I ) ; 






2 


V3-V0LTS( I NODE, 3 t I 3 : 


9 


L 


,\ 


V4=V0LTS( INOOE ,4 , I ) ; 


IV 


i 


-) 


OfcLTA_E( 1 tl )=V1-V^: 


1 1 


1 


I 


i)FL T*_E ( ? , I ) = V3-VA; 


I? 


1 


} 


DE LT A_E < 3 , 1 >=V2-V4 ; 


1 3 


1 


2 


0ELTA..E(4,I >=V3-Vl ; 


14 


1 


2 


OELTA_E (5,1 )=V^-V2; 
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SAE10: PROC OP T I ON S < MA r wi : 



STMT LEVE 



NEST 

2 
2 



1 
2 
2 
2 
3 
3 
3 
2 
2 
2 
2 
3 
3 
3 
2 
2 
2 
2 
2 
2 
2 
2 
2 
1 



LDU 



DELTA_E ( 6 1 1 >=V2-V1 ; 
END ; 

TRANSFORM KVA AND POWER FACTOR TO LOAD ADMITTANCE 
DO IY=1 TO 6; 

KVAZ=KVA< INODE, IV) ; 
IF KVAZ=0 
THEN DO; 

LOADY( INODE , IY,*)=C ; 
GO TO LDl; 
END ; 

PF=POWER_FACTGR( I NODE , IY) ; 
SINPF=SQRT(1-PF*PF) ; 

IF PFO 

THEN DO; 

SINPF=-SINPF ; 

pf=-pf; 

END ; 

E X = DE LTA_E ( I Y ? 1 ) ; 

EY=DELTA_E( IY,2 > ; 

EX=EX*EX; 

EY=EY*EY; 

ESUM=EX+EY; 

LODZ=KVAZ/ESUM; 

LOADY ( I NODE , I Y , 1 ) -LQDZ*PF ; 

LOADY( INODE, IY,2 > =L OD Z*S I NPF ; 

end; 

END; 



*/ 



/*- 



THIS SECTION REDUCES 8 NODE CIRCUIT TO 4 NODE EQUIVALENT 
DRIVING POINT NETWORK THROUGH STAP TO MESH TRANSFORMATIONS* 



-*/ 



1 

1 
1 

2 
2 

2 
2 
2 
2 
2 
? 
2 
3 
3 
2 



FOR ALGORITHM */ 



57 



DO I NODE =L A ST_NODE TO 2 BY -1; 
JNODE =FRGM_NODE ( I NUDE ) ; 
WORK=0; /* FILL WORKING ARRAY 

DO 1=1 TO 2; 

WORM It 1» I >=DPNY( I NODE, 4, I 

WORM 1.2 1 I ) -D PNY ( INODE ?5 ? I 

WORM 1 ?3, I )=DPNY( I NUDE, 3, I 

WORM 2,1,1) =D PNY ( I NODE ? 6 ? I 

WORM 2,2, I >=DPNY< I NODE, 2, I 

WORM 3, 1 , I ) =D PNY ( I NODE, 1,1 

W0RM4,f, I )=NEUTY( INODE, I ) ; 

PHASE =PHASEY< INODE, I ) ; 

DO J=l TO 3; 

WORM*, J, I )=PHASE ; 

end; 
end; 

START OF ALGORITHM FOR STAR - MESH TRANSFORMATION 
DO R=l TO 4: 
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STMT LEVEL NEST 



^9 

<r y . 
61 
6 2 

6 3 
64 
65 
66 
67 
68 
6<3 
7 ') 
71 
72 

7 3 
74 
7S 
76 
77 
7a 
79 
H 0 
ftl 
P2 



8 * 
84 
85 
86 

87 
R ft 
P9 

9 0 
91 
92 
93 
94 



2 

3 

3 

3 

2 

2 

2 

3 

3 

3 

4 

4 

4 

4 

4 

4 
4 
4 
4 
4 
4 
3 
2 



1 

2 
2 
2 
2 
2 
2 
2 
3 
3 
2 
1 



00 K=l TO 4; 

SUMX=SUMX+ WORK(K,R f l»; 
SUMY=SUMY+ WORK(K,R,2>; 

end; 

SUMSQ=SUMX*SUMX+ SUMY*SUMY; 

ri=r+i; 

DO S=2 TO 4; 
SM1-S-1 ; 
SMR1=S-R1 ; 
DO T=S TO 4; 

W0RK1X=W0RK(T-SM1 .Rtl ) 

WORKlY=WORK(T-SMl f Rt2) 

W0RK2X=W0RK(T,R,1) ; 

WORK2Y=WORK{T,R f 2) ; 

W0RK12X = W0PK1X**I0RK?X 

W0RK12Y-tfGRKlX*WQRK?Y 

PRODXM W0RK12X*SUMX «■ 

PRODY=( W0RK12Y-SUMX - 

TT=T-SMR1 ; 

WORK(SMl,TTfl»«PRODX* 
W0RK(SM1,TT,2)=PR0DY+ 

END ; 
END ; 

end; 

ADD DPN ADMITTANCES TO LOAD ADMITTANCES OF JNODE TO 
GET THE TOTAL DPN Y FOR JNODE 
DO 1=1 TO 2? 

DPN(1 WORKOtStI » ; 

0PN(2*« WORK( 2 » 6 ? I )\ 

DPN( 3)= WORK( 1 1 7 f I ) ; 

DPN(4)= WORM 1 f 5, f » ; 

DPN( 5 )= WORM 1 »6, I ) ; 

DPN( 61= WORK( 2 f 5t I> ; 

00 DP NY ( JNODE t IY» I ) =LQADY ( JNODE ? IY f I ) + DPN(IY); 

end; 
end; 
END; 



- WCRK1Y*W0RK2Y; 
+ W0RK1Y*W0PK2X; 
*ICRK12Y*SUMY)/SUMSQ; 
W0RK12X*SUMY»/SUMSQ; 

W0RK(SM1,TT,1»; 
WORK ( SMlt TT,2 » ; 



*/ 



'* "7c SECTION CALCULATES NODE VOLTAGES AT THE TO-NODE GIVEN 
TH^VO^AGES "THE FROM NODE f THE LINE ADMITTANCES BETWEEN 
tmp NHOFS AND THE DPN ADMITTANCES AT THE TO-NODE, THE 
SOLUTION I S OBTAINED BY SOLVING 8 SIMULTANEOUS REAL LINEAR 
EOS BY GAUSSIAN ELIMINATION WITH PARTIAL PIVOTING. THE EQS. 
ORIGINATE AS 4 COMPLEX NODE EQS ^=1^. *, 



95 
96 
97 
98 

99 



DO IN0DE=2 TO LAST_NODE ; 
JNODE=FROM_NODE ( INODE) ; 

DO 1 = 1 TO 2; /* RE AND I MAG PARTS */ 
Y1=DPNY( I NODE »1 * I ) ; 
Y2=DPNY( I NODE »2t 1 ) ; 
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PI SAfc'lO: PROC OP T I ONS { MA I N I ? 



S»MT L E V P L NEST 



ir i 

ir 2 

103 

i r a 
ins 



106 
K 7 

10 8 
10 q 

11^ 
111 

112 
113 
114 
115 
116 



117 
11« 

no 

12^ 
121 
122 
123 
124 



125 
126 
127 
128 
129 
130 
131 



2 
2 
? 
2 
2 
2 
2 
2 
2 
2 
? 



1 
1 
1 

2 
2 
2 
2 
2 



Y3=aPNY<INODE,3,I 
Y4=DPNY { I NODE ,4,1 
Y5=DPNY{ I NODE ,5,1 
Y6=DPNY< INOOE ,6,1 
YP=PHASEY< I NODE , I 
YN = NE UT Y{ I NOD E , I ) 

NODAL ADMITTANCE MATRIX */ 
Y{ 1,1 ,1 )=Y1 + Y4 + Y6 + YP; 
Y(2,2,I )=Y6+Y 3+ Y5+ YP ; 
Y<3,3,I )=Y4+Y2+Y5+YP ; 
Y(4,4,I )=Y1+Y3+Y2+YN; 
Y(1,2,I ),Y<2,1,I>=~Y6; 
Y( It 3 f I If Y(3tl t I ) =-Y4 ; 
Y{ 1,4,1 ) ,Y(4,1 ,1 ) =-Yl ; 
Y(2,3,I I ,Y( 3t 2 ,1 ) =-Y5 ; 
Y<4,2,I > ,Y<2,4,I >=-Y3; 
Y(4,3,I ),Y(3,4,I)=-Y2; 

end; 

/******** form current vfctor f ok eo ye* i */ 

Y PX=PHASEY( I NODE ,1 ) ; 
YPY=PHASFY( I NODE ,2) ; 
DO IROW=l TO 4; 

VOLTX=VOLTS< JNODf , I ROW, 1 ) ; 
VOLTY=VOLTS( J NODE ,IRJW,2) 5 
CUR( I ROW) =VOLTX*YPX-VOLTY*vpv ; 
CUR( I R0W+4)=V0LTX*YPY + V0LTY*YPX; 

end; 

/******** the 4 simultaneous complex linear eqs. are equivalent to 

b SIMULTANEOUS REAL LINEAR EOS, THESE ARE FORMED NEXT. 
THEY ARE YE=I < = > I YX -YY i I EX i - I I X J 

I YY YXI I EY| - I I Y| */ 

/******** FORM 8X8 MATRIX EQUIVALENT OF 4 CCMPLEX EQS, */ 

1 DO 1=1 TO 4; 

2 DO J=l TO 4; 

3 Z(I,J),Z(I+4,J+4)=Y(I,J,1) ; 

^ z(i+4 f J)«Y(if Jf2>; 

3 ti I f J+4)=-Y(I , J, 2) ; 

3 end; 
2 end; 



THIS ROUTINE DOES GAUSSIAN ELIMINATION WITH PARTIAL PIVOTING 
TO SOLVE THE 8X8 SYSTEM, 



-* r 



f if.if.if.lf.ifif.if if. LU DECOMPOSITION */ 

132 11 DO 1=1 TO 8; /* INITIALIZE POINTER FOR INTERCHANGING ROW*/ 

133 1 2 PT(I)=I; 
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ST LEVEL NEST 



1 

2 
2 
3 
3 
3 
4 
4 

3 
2 
2 
3 
3 
3 
3 
2 
2 
■a 

3 
3 
3 
4 
5 
5 
4 
3 
? 



1 

2 
2 
2 
3 
3 
2 
2 



1 

2 
2 
2 
3 
3 
2 
2 



DO K = l TO 7; 
MAX=0 ; 

DO I = K TO 8; /* SEARCH FOR MAX ELEMENT IN COL */ 
SIZE=ABS( Z(PT( I ) t K ) ) ; 
IF SIZE>MAX 
THEN DO; 
MAX=SIZE ; 
PI VOT_INDEX=I ; 
END ; 

end; 

IF PIVOT_INDtX-* = K 
THEN DO; 

SAVE_PT=PT( K) ; 
PT(K)=PT<PIVOT_INOEX) ; 
PTC PIVOT_INDEX)=SAVF_PT ; 
END; 

PIVOT=Z(PT<K> ,K) ; 
DO I=K + 1 TO 8 ; 
PTI=PT( I ) ; 

Z< PTI,K) ,MULT ^=Z(PTI ,K)/PT VOT; 
IF MULT-*=C 
THEN DO; 

DO J=K+1 TO B ; 

Z(PTI,J)=Z(PTI,J>-MULT*Z(PT(K> f J) ; 
END ; 
END; 

end; 
end; 

FORWARD EL IMJ NATION 
DU 1=1 TO 8; 
su M =c; 
PTI=PT< I ) ; 
DO J = l TO 1-1 ; 

SUM=SUM+ MULTIPLY<Z(PTI , J) f V(J), 16); 
END; 

V { I )=CUR( PTI ) -SUM ; 
END; 



/******** BACKWARDS SUBSTITUTION 
DO 1=8 TG 1 6Y -l; 
SUM=0; 
PTI=PT( I ) ; 
DO J=I+1 TO 8; 

SUM=SUM+ MULT I PLY <Z< PTI, J I •V(J),16) ; 
END ; 

V< I ) = < V ( I >-SUM)/Z<PTI ,1 ) ; 
END; 



*/ 



*/ 



179 1 1 IF TEST^=0 

18G 1 1 THEN GO TO NEWV; 



/******** TEST FOR MAGNITUDE OF CHANGE FROM PREVIOUS ITERATION */ 
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PI 


UE10: 


PROC OPTIONS! MAIN ) ; 


s r 


LFVhL 


NEST 






!. i 


1 


1 




DO 1=1 TO 2; 




1 


2 




IF 1 = 1 


i<n 


1 


2 




THEN IJ=0; 




1 


2 




ELSE IJ=4; 


1 35 


1 


2 




DO J=l TO 3; 


Lb6 


1 






vu=v< j + ij ) ; 


1?7 


1 






DI FFERENCE=ABS< ( VOLTS ( I NODE r J y I 


1 8 ^ 


1 


3 




IF DIFFERENCE>TOLERANCE 


L*9 




3 




THEN DO: 


1 r 5 '1 


j 


4 




TEST^l ; 


1 0 1 
1 / I 








GO TO NEWV; 


19? 




4 




end; 


193 








END ; 


194 


1 


? 




END ; 


195 


1 


1 


NEWV : 












DO 1=1 TO 4; 


196 




? 




VOLTS< I NODE , I , 1) -V( I ) ; 


197 




2 




VOLTS( INCOEt I »2>=V( 1+4) ; 


1<*R 




2 




end; 


199 




1 




end; 


2C0 


1 






IF TEST=1 


2M 


1 






T H FN GO TO ITERATE; 


? p 2 


1 






END PLSAEIO; 
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INSTRUCTOR'S GUIDE 



CASE TITLE : Radial Load Flow Problems in Distribution 
Circuit Analysis 

RELATED UNDERGRADUATE CURRIC ULUM A REAS : 

1. Engineering analysis and design methods (electri- 
cal engineering and computer science orientation) . 

2 . Electric circuit analysis (Node- voltage method , 
y-A transformations, practical application of 
power factor correction concept). 

3. Computer science 

a. Numerical- -iteration procedure, Gaussian elim- 
ination . complex numbers. 

b. Non-numerical- -traversing a tree data struc- 
ture in pre and end-order* 

4. Power systems analysis and control. 



SYNOPSIS I Calculation of the voltage drop from one point 
to the next point along a three-phase, low-voltage distri- 
bution circuit's transmission line can be done by approxi- 
mate means fairly accurately. However, in some cases, ex- 
act calculation of drops was found desirable by engineers 
at the Pennsylvania Power and Light Company. The problem's 
formulation introduces a large number of simultaneous non- 
linear algebraic equations which have complex co-eff icients- 
an obvious instance where the digital computer and numer- 
ical methods can be used to advantage. 

In this case, an electrical engineer and a scientific 
programmer work together to develop a computer program which 
depends on electric circuit equations and an iterative pro- 
cedure to solve them. The tree-like or radial nature of 
the typical distribution circuit required devising a way 
to sequence the computations properly. After the program 
was developed, a test was run on a simple circuit which 
could also be hand-calculated, and it was shown that com- 
puter results are quite different from the approximate re- 
sults which otherwise would have had to be accepted. 
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The simple concepts behind the iterative process and 
the computation sequence problem are easily shown to stu- 
dents without cluttering up the basic ideas with too much 
algebra. An interesting assignment in computer processing 
can thus be devised 

The case study also shows how locating capacitors at 
various points indicated by computation can improve volt- 
age regulation and lessen thermal losses in lines and 
transformers - -an illustration of the utilities 1 concern 
with avoiding waste of electric energy. 



QUESTIONS EQR THOUGHT Ad£ PIS QUSSJON 

After reading Part A: 

1. What are the primary reasons for controlling volt- 
age swings between no-load and full-load conditions 
on a distribution circuit? 



2. What part did a literature search play in the 
course of the project? 

3. A decision was made by the engineer on the extent 
to which he would explain the problem to the sci- 
entific programmer who was not an electrical engi- 
neer. Discuss this decision. 

4. Clearly explain the concept of power-factor cor- 
rection. Why might a large power system find it 
better to operate at 0.85 or 0.90 pf rather than 
nearer unity? 



After reading Part B: 



Given the problem below, solve for the two possible 
load voltages* What possible difficulty is sug- 
gested with respect to the way the more complex 
problem is solved in the case study? Discuss . 

i: 



+ o- 



125 V 



V L 



□ 

LOAD 

zj— 1 



2. 5 kw 
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Hint: The question raises the problem always re- 
lated to an iterative procedure- -that of conver- 
gence If the iteration procedure used in the case 
study is used here, the assumed starting value 
for will lead to convergence to one of the two 

possible values. 

Discuss the programming of the simpler DC problem 
below with the aim of making a general program, 
i.e., if a different set of cards as listed is 
used for input, the new circuit would be traversed 
and computed correctly. 




From Node To Node R(«) Load at To Node (w) 

0 1 .0010 2.5 

1 2 .0020 5.0 

2 3 .0005 10.0 

3 4 .0015 5.0 
2 5 .0020 6.0 
5 6 .0010 10.0 
5 7 ,0010 8.0 

Answers after 10 iterations (with 'flat' start) 
on IBM-360/65: 



V0 
VI 
V2 
V3 



1.0000 
.9412 
. 8290 
,8198 



V4 
V5 
V6 
V7 



.8105 

. 7655 
. 7522 
. 7549 



